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L  INTRODUCTION 


We  consider  multi-hop  packet  radio  topologies,  where  environmental  disturbances  dictate 
the  deployment  of  dynamic  assignments  for  the  transmission  channels  or  frequencies.  Then, 
dynamically  reconfigurated  cluster  topologies  result,  and  system  connectivity  is  ensured  via 
store-and-forward  operations  deployed  by  users  located  in  cluster  boundaries.  At  each  point  in 
time,  the  system  induces  a  topology  as  that  exhibited  in  Figure  1,  where  the  same  topology  has 
been  considered  in  [1].  In  each  cluster,  a  single  common  channel  is  used  for  transmissions  by 
the  ordinary  nodes  (or  users)  in  it.  The  clusterhead  in  each  cluster  is  one  of  the  ordinary  nodes 
(or  users)  and  its  function  is  to  monitor  the  multiple  access  transmission  algorithm  in  the  cluster, 
and  to  provide  the  required  feedbacks.  The  gateway  nodes  together  with  the  clusterheads  form 
the  backbone  network  in  the  system,  and  ensure  intracluster  connectivity.  The  dynamics  of  the 
cluster  formations  and  the  structure  of  the  backbone  network  are  directly  related  to  the  opera¬ 
tional  and  performance  characteristics  of  the  deployed  multiple-access  algorithms. 

In  the  topologies  considered  here,  the  performance  criteria  include  reliability,  connectivity, 
economy,  bandwidth  efficiency,  effective  control  of  transmission  delays,  and  robustness  in  the 
presence  of  channel  errors.  For  reliability,  economy,  and  bandwidth  efficiency,  the  deployed 
multiple  access  algorithms  should  require  very  little  knowledge  about  the  system  state  by  the 
users,  since  such  knowledge  can  be  obtained  only  at  the  expense  of  available  bandwidth  and  its 
use  requires  complex  and  expensive  hardware,  and  since  the  mobile  users  are  generally  unable  to 
maintain  overall  system  state  information  at  all  times.  Thus,  for  the  present  topologies,  we  only 
consider  Limited  Sensing  Random  Access  Algorithms  (LSRAAs)  for  transmission,  (as  those  in 
(3]  and  [4]).  Those  algorithms  only  require  that  each  user  know  the  algorithmic  operations,  can 
detect  the  transmission  frequency  of  the  cluster  he  is  located  in,  and  can  sense  the  channel  feed¬ 
back  from  the  time  he  generates  a  packet  to  the  time  that  this  packet  is  successfully  transmitted; 
they  do  not  require  additional  knowledge  about  the  system  state  or  topological  information. 

As  the  topology  in  Figure  1  evolves  dynamically,  and  for  better  system  connectivity,  neigh¬ 
boring  clusters  may  overlap  (see  [1],  for  example).  The  ordinary  nodes  locate-  the  overlap¬ 
ping  regions  are  then  exposed  to  transmissions  and  feedbacks  from  more  than  one  clusters;  a 
generally  time-varying  phenomenon  due  to  the  mobility  of  the  ordinary  nodes,  which  can  be 
exploited  for  the  improvement  in  performance  of  the  overall  system,  (as  done  in  [2],  for  the 
mobile  telephone  system).  Consider,  for  example,  clusters  1  and  2  in  Figure  1,  and  let  us  call  the 
ordinary  nodes  in  their  overlapping  region,  marginal  users;  let  us  call  the  ordinary  nodes  in  clus¬ 
ter  i,  i=l,2,  which  are  not  located  in  the  overlapping  region,  local  users  for  cluster  i.  The  local 
users  in  cluster  i  communicate  via  the  LSRAA  deployed  by  the  cluster,  called  LSRAAi,  i=l,2. 
Due  to  the  double  exposure,  the  marginal  users  have  a  choice:  For  communication  with  each 
other  and  system  connectivity,  they  can  join  either  one  of  the  LSRAAI  and  LSRAA2  algo¬ 
rithmic  systems.  This  choice  can  be  implemented  either  dynamically,  dictated  by  the  two  feed¬ 
back  sequences  (from  the  two  clusters)  that  the  marginal  users  observe,  or  statically;  the  specif¬ 
ics  of  its  implementation  will  be  called,  interconnection  algorithm.  A  static  interconnection  algo¬ 
rithm  is  represented  by  a  priori  assigned  probabilities.  Then,  upon  generation  of  a  new  packet, 
each  marginal  user  joins  the  LSRAAi  with  probability  p;,  i=l,2,  where  Pi  +  P2  =  1,  and  remains 
there  until  his  packet  is  successfully  transmitted.  The  above  static  policy  is  simple,  but  requires 
that  the  marginal  users  know  the  a  priori  assigned  probabilities  at  all  times.  In  the  dynamically 
changing  topologies  considered  here,  those  probabilities  should  change  dynamically  and  their 
values  should  remain  knowr  to  all  users.  But  this  implies  kuo-.vledge  of  the  system  dynamics  at 
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all  times,  which  is  either  very  hard  to  obtain  or  requires  a  tremendous  increase  in  feedback  infor¬ 
mation,  (and  thus  tremendous  increase  in  the  bandwidth  of  the  feedback  channels).  Therefore, 
only  dynamic  interconnection  algorithms  are  appropriate  for  the  mobile  multi-hop  packet  radio 
system.  Such  algorithms  may  in  addition  give  an  advantage  to  the  marginal  users,  in  terms  of 
delays.  This  may  be  desirable  in  systems  where  the  marginal  users  transmit  either  priority  mes¬ 
sages  or  control  data,  (including  topological  information  for  reconfiguration  of  the  cluster  struc¬ 
ture,  etc.). 

In  this  paper,  we  consider  two  overlapping  clusters,  each  deploying  a  LSRAA  for  its  local 
users.  Then,  we  propose  a  dynamic  interconnection  algorithm,  via  which  marginal  users  join 
either  one  of  the  local  per  cluster  LSRAA.  For  each  of  the  two  local  LSRAAs,  we  adopt  the  lim¬ 
ited  sensing  form  of  the  algorithm  in  [4].  The  latter  algorithm  has  simple  operational  properties, 
in  the  presence  of  the  limit  Poisson  user  model  it  attains  throughput  0.43,  it  has  superior  resis¬ 
tance  to  feedback  errors,  and  operates  with  binary-collision  versus  noncollision-feedback. 

The  organization  of  the  paper  is  as  follows:  In  Section  II,  we  present  the  system  model.  In 
Section  m,  we  describe  the  local  LSRAAs  and  the  algorithm  that  the  marginal  users  deploy.  In 
Section  IV,  we  present  the  algorithmic  analysis  of  the  two-cluster  system.  In  Section  V,  we 
include  comments  and  conclusions. 


II.  SYSTEM  MODEL 

We  consider  the  two-cluster  packet-radio  system  in  Figure  1.  We  assume  that  in  each  of 
the  two  clusters,  some  synchronous  LSRAA  is  deployed.  In  particular:  (1)  Time  is  divided  in 
slots  of  length  equal  to  the  duration  of  a  packet,  and  the  starting  instants  of  the  slots  are  identical 
in  both  clusters.  (2)  In  each  cluster,  the  clusterhead  broadcasts  a  feedback  per  slot,  which 
corresponds  to  the  outcomes  induced  by  the  local  LSRAA.  This  feedback  is  either  ternary.  Col¬ 
lision  (C)  versus  Success  (S)  versus  Emptiness  (E),  or  binary.  Collision  (C)  versus  Non-Collision 
(NC).  (3)  In  each  cluster,  each  local  user  is  required  to  monitor  the  feedback  'm  the  local 
clusterhead  continuously,  from  the  time  he  generates  a  new  packet  to  the  time  tha:  his  packet  is 
successfully  transmitted.  We  assume  that  no  propagation  delays  and  no  forward  or  feedback 
channel  errors  exist  in  the  system. 

We  assume  that  each  marginal  user  receives  the  feedbacks  from  both  the  local  clusterheads 
correctly  and  without  propagation  delays.  At  the  time  when  a  marginal  user  generates  a  new 
packet,  he  starts  monitoring  the  feedbacks  from  both  clusterheads  continuously,  until  he  decides 
to  join  the  operations  of  one  of  the  two  LSRAAs,  for  the  transmission  of  his  packet.  Upon  this 
decision,  he  maintains  the  continuous  monitoring  of  only  those  feedbacks  that  correspond  to  the 
LSRAA  he  chose,  until  his  packet  is  successfully  transmitted. 

It  is  assumed  that  the  mobility  of  the  users  in  the  system  is  low  enough,  so  that  each  user 
remains  within  the  same  geographical  region  (local  or  marginal)  from  the  time  he  generates  a 
packet  to  the  time  that  this  packet  is  successfully  transmitted.  If  the  local  LSRAAs  in  each  clus¬ 
ter  have  good  delay  characteristics,  then  this  time  period  may  be  relatively  small,  with  high  pro¬ 
bability. 

The  local  user  populations  in  each  cluster  and  the  marginal  user  population  are  all  modelled 
as  limit  Poisson.  That  is,  it  is  assumed  that  the  local  traffic  generated  in  cluster  i,  i=l,  2,  is  a 
Poisson  process  with  intensity  ,  i=  1 ,2,  and  that  the  traffic  generated  by  the  marginal  users  is 


another  Poisson  process  with  intensity  X3 .  As  found  in  [5],  for  a  large  class  of  LSRAAs,  the 
limit  Poisson  user  model  provides  a  lower  bound  in  performance,  within  the  class  of  identical 
and  independent  users  whose  packet  generating  process  is  i.i.d. 

EL  THE  ALGORITHMS 

We  assume  that  the  two  LSRAAs  in  the  system  are  identical.  Each  LSRAA  is  the  window 
algorithm  in  [4],  which  induces  throughput  0.43  and  operates  with  binary  C  versus  NC  feedback. 
This  algorithm  has  simple  properties  and  operations,  while  the  algorithm  in  [3]  is  more  complex. 
In  addition,  the  algorithm  in  [4]  has  superior  resistance  to  feedback  errors. 

Upon  generation  of  a  new  packet,  a  marginal  user  imagines  himself  belonging  to  the  sys¬ 
tems  of  both  the  LSRAAs  and  follows  their  algorithmic  steps,  until  the  first  time  that  he  enters  a 
collision  resolution  event  in  one  of  them.  Then,  he  remains  with  the  latter  LSRAA  system,  until 
his  packet  is  successfully  transmitted. 

In  this  paper,  we  consider  the  case  where  each  local  LSRAA  is  the  limited  sensing  version 
of  the  two-cell  algorithm  in  [4].  The  reasons  for  this  choice  are  several:  This  algorithm  attains 
high  enough  throughput  and  induces  low  delays,  while  at  the  same  time  has  simple  operational 
and  analytical  properties,  and  is  very  resistant  to  feedback  errors.  For  completeness,  we  describe 
the  algorithm  here. 

Let  time  be  measured  in  slot  units,  where  slot  t  occupies  the  time  interval  [t,  t+1).  Let  xt  (j) 
denote  the  feedback  that  corresponds  to  slot  t,  for  cluster  j;  j=l,2,  where  xt(j)=C  and  xt(j)=NC 
represent  collision  and  noncollision  slot  t  in  cluster  j,  respectively.  The  local  LSRAA  in  cluster  j 
is  implemented  independently  by  each  user  in  the  system,  and  utilizes  a  window  of  length  A  Let 
some  local  in  cluster  j  user  generate  a  new  packet  within  the  time  interval  [tj ,  ti+1).  Then,  he 
immediately  starts  observing  the  feedback  sequence  ,  beginning  with  the  feedback 

xt(  (j)-  Let  us  define  the  sequence  {tjQ}^,  as  follows:  t2<J)  *s  the  first  time  after  :t,  such  that 
xt2(j)(i)  =  xt2<j>— 1  =  NC.  Then,  as  will  be  explained  below,  t20)  corresponds  to  the  _.-.^ing  slot  of  a 
Collision  Resolution  Interval  (CRI)  in  cluster  j,  and  from  t2(j)+l  on,  the  user  can  identify  the 
ending  slots  of  CRIs  induced  by  the  algorithm  in  cluster  j.  Each  t;(j)  corresponds  to  the  ending 
slot  of  some  CRI  in  cluster  j,  and  ti+j  (j)  is  the  first  after  t;(j)  such  slot.  At  tj(j)>  the  user  updates 
his  arrival  instant,  as  follows:  tjl)  =  tj  +  (i-2)A;  we  call  the  sequence  {t^}^,  updates.  Let  tk(j) 
be  such  that:  tk(j)e(ti(j)}i>2>  tr<ti(j)-l-A  ;  Vi<k-1,  and  t^k)>tk(j)-l-A.  Then,  in  slot  tk(j)+l, 
the  user  enters  a  CRI  within  the  LSRAA  of  cluster  j,  and  transmits  his  packet  successfully  during 
its  process.  He  stops  observing  the  feedback  sequence  { xt(J) }  at  the  point  when  his  packet  is 
sucessfully  transmitted.  If  the  user  is  instead  marginal,  then  he  observes  both  feedback 
sequences  (xt(l)}t^tl  and  {xt(2)}ctl,  and  follows  the  evolution  of  both  the  time  sequences 
{t(l))i>2  and  {h(2)}&2-  If  tk(l)<tk(2),  then  in  slot  tk(l)+l  he  enters  a  CRI  within  the  LSRAA  of 
cluster  1,  and  transmits  his  packet  successfully  during  its  process.  If  tk(l)>tk(2),  instead,  then  he 
joins  a  CRI  in  cluster  2,  in  slot  tk(2)+l.  If  tk(l)=tk(2),  then  he  selects  one  of  the  local  LSRAAs 
with  probability  0.5.  The  above,  describe  the  first  entry  rules,  for  the  local  and  the  marginal 
users;  that  is,  how  and  when  each  newly  generated  packet  first  starts  participating  in  some  CRI, 
for  its  successful  transmission.  From  the  first  entry  rule  that  the  marginal  users  utilize,  it  is  clear 
that  they  have  an  advantage  over  the  local  for  clusters  1  and  2  users.  In  particular,  their  waiting 
time  until  they  first  enter  some  CRI  is  generally  smaller  than  that  of  the  local  users;  thus,  their 


58 


overall  delays  are  generally  smaller  than  those  of  the  local  users  as  well. 

Consider  the  algorithm  in  cluster  j,  and  let  it  start  operating  at  time  zero.  Then,  slot  1  is 
empty.  In  slot  2,  the  arrivals  in  [0,1)  are  transmitted,  and  a  CRI  begins.  If  the  number  of  arrivals 
in  [0,1)  is  less  than  two,  then  X2(j)=NC,  the  CRI  lasts  one  slot,  and  a  new  CRI  begins  with  slot  3. 
If  the  number  of  arrivals  in  [0,1)  is  at  least  two,  then  X2(j)=C,  instead,  and  the  CRI  lasts  as  long 
as  it  takes  to  resolve  the  collision  in  slot  2;  its  end  is  identifiable  by  all  the  users  in  the  system, 
(as  will  be  seen  below).  In  general,  let  T  be  a  slot  that  corresponds  to  the  end  of  some  CRI. 
Then,  in  slot  T+l,  all  the  users  with  current  updates  in  (T-A-1,T-1]  transmit.  If  Xx+i(j)=NC, 
then  the  CRI  which  started  with  slot  T+l  lasts  one  slot,  and  a  new  CRI  starts  with  slot  T+2.  If 
xT+i  (j)=C,  instead,  then  a  collision  occurs,  whose  resolution  starts  with  slot  T+2.  No  arrivals 
that  did  not  participate  in  the  collision  at  T+l  are  transmitted,  until  the  latter  is  resolved.  During 
the  collision  resolution,  each  involved  user  acts  independently,  via  the  utilization  of  a  counter 
whose  value  at  time  t  is  denoted  rt.  The  counter  values  can  be  either  1  or  2,  and  they  are  updated 
and  utilized  according  to  the  rules  below. 

1.  The  user  tranmits  in  slot  t,  if  and  only  if  rt=l.  A  packet  is  successfully  transmitted  in  t,  if 
and  only  if  rt=l  and  xt=NC. 

2.  The  counter  values  transition  in  time  as  follows: 

(a)  If  xt_i=NC  and  rt_!=2,  then  rt=l 

(b)  If  xt_!  =C  and  rt_i  =2,  then  rt=2 

(c)  If  xt_i  =C  and  rt_j = 1 ,  then 

{1  ,  with  probability  0.5 
2  ,  with  probability  0.5 

A  CRI  which  starts  with  a  collision,  ends  when  it  becomes  known  to  all  users  that  the  initially 
collided  packets  have  been  successfully  transmitted.  From  the  operations  exhibited  above,  it  is 
not  hard  to  see  that  such  a  CRI  ends  the  first  time  (after  its  beginning)  that  two  consecutive  NC 
slots  occur. 


IV.  ALGORITHMIC  ANALYSIS 

For  convenience  in  notation,  we  will  refer  to  the  local  users  in  cluster  1,  the  local  users  in 
cluster  2,  and  the  marginal  users,  as  subsystem  1,  subsystem  2,  and  subsystem  3,  respectively.  In 
the  algorithmic  analysis,  we  will  adopt  the  limit  Poisson  user  model  (infinitely  many  indepen¬ 
dent  Bernoulli  users)  for  each  of  the  three  subsystems.  In  particular,  we  will  assume  that  the 
three  subsystem  traffics  are  mutually  independent,  and  that  the  user  traffic  in  subsystem  j, 
j= 1,2,3,  is  limit  Poisson  with  intensity  Xj. 

Consider  either  one  of  the  local  LSRAAs,  whose  operations  are  described  in  section  III. 
Consider  some  CRI  within  the  system  of  the  LSRAA,  which  starts  with  transmissions  from 
packet  arrivals  in  an  arrival  interval  of  cumulative  length,  u.  We  will  call  u,  the  "length  of  the 
examined  interval."  Let  us  then  define: 

E(/ 1  u } :  Given  length  of  the  examined  interval  equal  to  u,  the  expected  number  of  slots 

needed  for  its  resolution;  that  is,  for  the  successful  transmission  of  all  the  arrivals 
in  the  examined  interval. 


a  a  a  -•  j ■«>>.> )>-> j.'/' *'«»  * r.'  •  V  V 


^%k-n: 


Given  that  n  packets  have  counter  values  equal  to  1,  and  k-n  packets  have  counter 
values  equal  to  2,  the  expected  number  of  slots  needed  by  the  local  LSRAA,  for 
the  successful  transmission  of  all  the  k  packets. 

Then,  directly  from  the  results  in  [4],  we  have: 


Lo.o  =  Li.o  =  1  »  Lo.i  =  1+Li.o  ;  £1 

0  <  L^q  <  jk2  +  j  k-2  ;  k>l 


(1) 


If  the  arrivals  in  the  examined  interval  are  controlled  by  a  limit  Poisson  process  whose 
intensity  is  X,  then: 

E{/lu}  *  Hx{/lu}  =  £  !*.<)  e-Xu  (2) 

fc=0  K- 

As  found  in  [4],  in  the  presence  of  the  limit  Poisson  user  model,  each  local  LSRAA  has 
throughput  0.43,  attained  for  window  size  A*  =  2.33.  Thus,  in  view  of  the  system  considered  in 
this  paper,  the  Poisson  intensities  Xj,  j=l,2,3,  must  satisfy  the  following  necessary  conditions,  for 
overall  system  stability: 

X!  <  0.43,  and  X2  <0.43,  and  Xt  +  X2  +  X3  <  2(0.43)  =  0.86  (3) 

The  necessary  conditions  in  (3)  determine  a  {Xj,j=l,2,3}  hyperplane,  which  contains  the 
{Xj,  j=l,2,3)  region  that  determines  the  system  throughput.  Tight  bounds  on  the  {Xj,  j=l,2,3) 
space  which  provides  the  system  throughput,  will  be  attained  via  the  system  stability  analysis  in 
Section  IV.  1  below. 


IV.l.  System  Stability 


We  consider  the  evolution  of  the  algorithms  in  the  two-cluster  system,  and  we  assume  that 
the  system  starts  operating  at  time  zero.  Let  us  consider  the  sequence  in  time  of  the  CRIs 
induced  by  the  two  LSRAAs  in  the  system.  Let  the  sequence  {TnJn^o  be  such  that:  (1)  For  each 
n,  Tn  corresponds  to  the  starting  point  of  a  slot  which  is  the  beginning  of  some  CRI.  We  note 
that  at  Tn,  two  CRIs  may  simultaneously  begin;  one  for  each  of  the  two  LSRAAs  in  the  system. 
(2)  Tn  is  the  first  after  Tn_i  time  instant  which  corresponds  to  the  beginning  of  some  CRI.  (3) 
To=2,  and  at  Tq  two  CRIs  begin;  one  for  each  of  the  two  LSRAAs  in  the  system. 


Let  {TnS)}n>o  be  the  subsequence  of  sequence  { Tn } „>o ,  which  consists  of  those  time 
instants  when  two  CRIs  begin  simultaneously;  one  for  each  of  the  two  LSRAAs  in  the  system. 
Clearly,  Tq  )=T0=2.  Let  D®s+1,  j=l,2,3,  denote  the  total  length  of  the  unresolved  arrival  inter¬ 
vals  in  subsystem  j,  at  the  time  instant  T(ns) .  D^s  is  then  called  "the  lag  of  subsystem  j  at  time 
Tls)."  From  the  algorithmic  operations  in  the  system,  we  conclude:  (1)  Dj^s>l  and  the  values  of 
D^j  are  denumerable  for  all  n  and  j.  (2)  D(jPs=I,  j=l  ,2,3.  (3)  At  time  t£s\  the  LSRAA  in  clus¬ 
ter  k,  k=l,2,  examines  two  arrival  intervals:  one  from  subsystem  k  which  has  length 
min(D^,  A)  and  contains  arrivals  generated  by  a  Poisson  process  with  intensity  Xk,  and  one 
from  subsystem  3  which  has  length  min  (D^,  A)  and  contains  arrivals  generated  by  a  Poisson 
process  with  intensity  0.5  X3.  (4)  The  triple  (D^s,  j= 1,2,3)  describes  the  state  of  the  system  at 

time  T„s),  and  the  sequence  {Snl^o  ^  {D^s,  j=l,2,3)n2o  is  a  three-dimensional  irreducible  and 


aperiodic  Markov  Chain. 

The  stability  of  the  system  is  represented  by  the  ergodicity  of  the  three-dimensional  Mar¬ 
kov  Chain  { Sn  Jn^o-  In  addition,  at  time  Tjf\  the  backlogs  of  each  of  the  two  LSRAAS  and  of 
the  overall  system  are  all  respresented  by  the  three  lags,  D^s,  j= 1,2,3.  Given  D^s  =  d§\ 

j=l,2,3,  the  expected  system  backlog  at  time  is  £  Xjdjjp,  where  V(d„)  = 

3 

V({d^)i^j<3) ^  £  \j  djp  is  a  Lyapunov  function  of  the  three  system  lags  {djp)i<js3.  Let  C 
j=i 

denote  the  state  space  of  the  Markov  Chain  {Sn  J^,  and  let  us  then  define  the  operator 
AV(d)  =  AV((d(j)  }1Sj<3),  called  a  generalized  drift,  as  follows: 

AV(d)  =  AV({d«)i<j<3)  =  H{V({D©u}1<j<3)-VaD®  }lsjs3)ID®s=dO)j=i,2,3} 

=  E{£xjD<&1.s  ID»=d®,j=l,2,3}  -  £xjd(j);d={d(j)}iSj<3  eC  (4) 

j=  i  i=i 

The  generalized  drift  in  (4)  can  be  used  to  establish  necessary  and  sufficient  conditions  for  the 
ergodicity  of  the  Markov  Chain  {Sn  }„><).  To  see  that,  let  us  define: 

P4  5_P(d„+i  =eldn=d)  (5) 


v;-,V(e)  . 


Fd.v(  Z^d-zr1  ZV(d)-IPieZ 


;  z  e[0, 1),  d  eC  (6) 

;  where  F^v  (z)  is  the  generalized  Kaplan  function  on  the  Lyapunov  function  V,  [11]. 
From  the  results  in  [7],  [8],  [9],  and  [10],  we  can  then  express  the  following  Proposidon. 

Proposition 


If  it  can  be  established  that: 

(i)  I  AV(d)  I  <°°,  for  all  d  in  C. 

(ii)  Theio  exist  finite  subset  H2  of  C  and  some  nonnegative  finite  constant  B,  such  that, 

F^v  (z)  ^  -B,  for  all  z  in[0, 1),  and  all  d  in  C-H2  (8) 

Then, 

(A)  If  there  exists  finite  subset  Hi  of  C,  and  some  e'>0,  such  that, 

AV(d)<-e',  for  all  d  in  C-Hj  (9) 

then,  the  Markov  Chain  {Sj,}^  is  ergodic. 

(B)  If  H  is  a  proper  subset  of  C,  such  that, 

_  inf  V(d)  >  sup  V(d) 


AV(d)  >  0,  for  all  d  in  C-H 
then,  the  Markov  Chain  { Sn  is  nonergodic.  □ 


(ID 


Remark.  Condition  (8)  is  trivially  satisfied  for  Markov  Chains  which  are  downward  uniformly 
bounded  with_  respect  to  the  Lyapunov  function  V;  that  is,  for  Chains  such  that, 
P(V(Dn+i)<V(dn)-x)=0;  for  all  x>a,  for  some  positive  and  finite  constant  a,  and  for  all  d,,  in  C. 
In  the  problem  studied  in  this  paper,  when  ^3=0,  two  decoupled  systems  arise.  Each  of  those 
systems  corresponds  then  to  a  local  LSRAA,  whose  lags  at  collision  resolution  points  form  a 
downward  uniformly  bounded  Markov  Chain;  the  Lyapunov  function  V(d)  reduces  then  to  d^\ 
for  the  system  in  cluster  j,  j=l,2,  and  the  constant  a  equals  then  [  Aj  =  2,  (see  [4]). 

We  now  express  a  lemma,  whose  proof  is  in  the  Appendix. 

Lemma  1 


Consider  the  Markov  Chain  {Sninso  =  {D^,  j=l,2,3)n2o  and  its  Lyapunov  function 

V(d„)  =  £  Xj  d$\  and  let  C  be  the  state  space  of  {Sn)n>o.  Let  AV(d)  and  be  as  in  (4) 

l<j£3 

and  (6),  respectively.  Then, 

_  3 

(i)  I A V(d)  I  <  oo.  for  all  d  in  C  and  all{Xj }:  £  ^1 

j=i 

(ii)  There  exist  positive  finite  constant  B  and  some  finite  subset  H2  of  C.  such  that: 


(12) 


F^vU)  -  ~B,  for  all  zin[0,l),  alldin  C-H2,andall  {X.j}:  £Xj  <  1  □ 

j=i 


(13) 


Due  to  Lemma  1,  and  in  view  of  the  Proposition,  the  search  for  the  ergodicity  versus  the 
nonergodicity  of  the  Markov  Chain  {Sn  reduces  to  the  identification  of  appropriate  proper¬ 
ties,  for  the  satisfaction  of  conditions  (A)  and  (B)  in  the  Proposition,  respectively. 

Consider  T^s),  and  let  the  lags  {D^s}i<j<3  be  sufficiently  long,  so  that  in  the  time  interval 
[T^s),  Tn+i],  each  CRI  from  each  of  the  two  LSRAAs  in  the  system,  resolves  an  arrival  interval 
of  length  A,  (as  proven  in  the  Appendix,  the  lengths  of  the  lags  {D^s }  i<j<3  that  satisfy  this  con¬ 
ation  are  finite  with  probability  close  to  one).  Given  system  Poisson  rates  { ^-j }  i<j<3 ,  and  for 
Ex  ( / 1  u }  as  in  (2),  let  us  then  define: 


i=l,2;Ni«\j},A): 


AV(dA,  {Xj }): 


Given  { Xj } ,  the  number  of  CRIs  generated  by  the  algorithm 
in  cluster  i,  in  [T^s),  T^t],  when  each  CRI  resolves  arrival 
intervals  of  length  A,  from  both  its  local  and  the  marginal 
users.  The  first  such  CRI  starts  at  T^\  and  the  last  such 
CRI  ends  at  T&. 

The  generalized  drift  in  (4),  for  given  Poisson  rates  {X,}, 
when,  in  jT^s),  T£.\],  each  CRI  generated  by  either  algo¬ 
rithm  in  the  system,  resolves  arrival  intervals  of  length  A. 


i=l,2;  A(i)V(dA,  {Xj})^  Ex^+Xj  12  {/IA)-A  + 


E  { Ni  ( { Xj } ,  A) } — 1 


Ex^Xj  { / 1 A }  — A 


(14) 


'« '* .  A  . 


.  A  AV. 
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Then,  via  (4),  we  easily  derive  the  following  expression: 

X  X 

AV(dA,  aj})  =  (Xl+-j-)A(1)V(dA,  {Aj})  +  (X2+-y-)A(2)V(dA,{Xj})+ 


X3  A 

~T~ 


E{N1({X]},A)}+E{N2({X3},A)}-2 


;  where, 

Ex1+X,i2{/IA}+[E{N1({Xj},A)}-l]EXl+x,{/IA}= 


(15) 


=  Ex2+x,i2{/!A}  +  [E{N2({Xj},A))-l]EX2+X3{/IA}  (16) 

We  note  that  AwV(dA,  {Xj })  in  (14)  is  a  generalized  drift  for  the  algorithm  in  cluster  i.  We  now 
state  a  theorem  whose  proof  is  in  the  Appendix. 


Theorem  1 

(i)  Let  there  exist  some  e>0,  such  that  the  two  conditions  below  are  satisfied: 

A(1)V(dA,{Xj})  <-e 

A(2)V(dA,  {Xj})  <-e  (17) 

Then,  the  Markov  Chain  (Sn)n>o  is  ergodic  at  the  Poisson  rates  {Xj  h<j<3  . 

(ii)  Let  at  least  one  of  the  two  generalized  drifts,  A(l)V(dA,  { Xj  }),i=l,2,  be  nonnegative. 
Then,  the  Markov  Chain  {Sn  )n>o  is  nonergodic  at  the  Poisson  rates  { Xj }  i^j<3  •  □ 

Remark:  The  stability  analysis  remains  unchanged,  when  the  two  random  access  algo¬ 
rithms,  in  clusters  1  and  2,  are  full  feedback  sensing,  instead.  The  full  versus  limited  feed¬ 
back  sensing  choice  affects  only  the  transmission  delays. 

The  conditions  (17)  in  the  Theorem  define  lower  bounds  on  the  (Xj)  regions  for  which  the 
Markov  Chain  {Sn)n>o  is  ergodic.  Since  the  constant  E  in  (17)  can  be  arbitrarily  small,  and  in 
combination  with  the  statement  in  part  (ii)  of  the  Theorem,  we  conclude  that  the  ergodicity 
bounds  determined  by  (17)  are  tight. 

In  the  Appendix,  and  in  Section  VI.2  in  particular,  we  describe  the  methodology  regarding 
the  computation  of  the  {Xj)  values  for  which  the  Chain  {Sn)n>o  is  ergodic,  as  well  as  the  per¬ 
tinent  quantities  and  their  recursions.  In  Table  1,  we  include  the  maximum  rate  X2  accepted  by 
the  system,  for  varying  Xi  and  X3  values,  as  well  as  the  maximum  accepted  symmetric  rate 
X*(X3)  =  X*  =X*  =  X2,  for  varying  X3  values.  We  also  include  the  maximum  rate  X3  accepted 
by  the  system,  when  Xi  =  X2  =  0.  We  note  that  due  to  symmetries  appearing  in  the  system  when 
XtandX2  are  interchanged,  for  each  fixed  rate  X3,  we  only  need  to  compute  X2  for  X[  values  in 
the  interal  [0,X*(X3)).  In  Figure  2,  we  plot  the  boundaries  of  the  (Xj,  X2)  acceptable  regions, 
parametrized  by  various  X3  values.  Those  boundaries  are  clearly  symmetric  around  the  45° 
straight  line.  In  Figure  3,  we  plot  X3  against  X  =  X3  =  X2,  for  symmetric  two-cluster  systems. 

From  Table  1  and  Figures  2  and  3,  we  observe  that  whenever  X3  is  strictly  positive,  then  the 
sum  X*+X2+X-  is  strictly  less  than  0.86;  that  is,  strictly  less  than  twice  the  throughput  of  each 
local  LSRAA.  This  is  so,  because  when  X3>0,  some  percentage  of  the  traffic  generated  by  the 
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marginal  users  is  always  assigned  for  transmission  to  an  even  over'o^ded  by  local  traffic 
LSRAA.  Thus,  if  the  rate  of  the  local  traffic  is  close  to  the  algorithmic  throughput,  then  the 
presence  of  marginal  users  results  in  instability  of  the  LSRAA. 


IV .2.  System  Delays 

The  two-cluster  algorithmic  system  induces  regenerative  points  within  its  stability  region. 
Thus,  in  essence,  the  methodology  in  [6]  applies  for  the  computation  of  the  expected  delays,  in 
each  of  the  three  subsystems.  To  see  that,  consider  the  sequence  {Sn}n^o  —  {D^s,  j=l,2,3}n20  of 
lag  vectors  induced  by  the  algorithmic  system,  at  the  time  instants  {T^}„>o.  As  established 
before,  the  lag  vector  So  is  the  unity  vector.  Let  us  define  the  sequence  {Ri)i>i  as  follows: 
Rt  =  Tqs)=2,  and  each  R,  corresponds  to  some  T^s)  instant  such  that  the  lag  vector  Sn  is  the  unity 
vector,  Ri+l  is  the  first  after  Ri  such  instant.  Let  QP,  i>l,  j=  1,2,3,  denote  the  number  of 
transmitted  packets  from  subsystem  j,  in  the  time  interval  (  0,R;].  Then,  GP^QPi  -QP 
denotes  the  number  of  transmitted  packets  from  subsystem],  in  the  time  interval  (Rj,  Ri+1],  and 
for  memoryless  traffics,  such  as  the  Poisson,  the  sequences  {GP  J^i,  j=l,2,3,  are  sequences  of 
i.i.d.  random  variables;  thus,  the  sequences  {QP)i>i,  j=l,2,3,  are  then  renewal  processes.  In 
addition,  if  DP  denotes  the  delay  experienced  by  the  n-th  transmitted  packet  arrival  from  sub¬ 
system  j,  then,  the  delay  process  {Dp  }„>i ,  j=1^2,3,  induced  by  the  algorithmic  system  is  regen¬ 
erative  with  respect  to  die  process  {QP)i>i,  and  the  process  (GP  )i>j  is  nonperiodic  for  every  j, 
since  P(Gp=I)>0,  j=  1,2,3.  Let  D(j)  denote  the  expected  steady-state  delay  experienced  by  a 
packet  in  subsystem  j,  j=l,2,3,  and  let  us  define: 

Z(i)  ^  E{GP)  ,  j=l ,2,3 
a  G? 

DPl  .  j=l,2,3  (18) 

i=l 

H  ^  E{R2  -  Rj ) 

Then,  from  the  regenerative  arguments  in  [6],  we  conclude: 

D(i)=W(i)  [Z<i>  rl.j=l,2,3  (19) 

where,  Z®  =  Xj  H ,  j=  1,2,3  (20) 

The  computation  of  bounds  on  the  expected  delays  depend  on  the  computation  of  tight 
upper  and  lower  bounds  on  the  quantities  in  (18).  The  pertinent  quantities  towards  that  direction 
and  their  recursions  are  included  in  the  Appendix.  We  used  them  together  with  the  methodology 
in  [6]  to  compute  expected  per  packet  delays.  We  present  some  of  our  results  in  Figures  4,  5,  6, 
and  7.  In  Figures  4  and  5,  we  plot  expected  per  packet  delays,  for  two  symmetric  systems 
(X]  =  X2),  as  functions  of  the  Poisson  intensity  X3.  In  Figures  6  and  7,  we  plot  expected  delays 
against  X3,  for  two  asymmetric  cases.  From  the  four  figures,  we  observe  the  advantage  of  the 
marginal  users,  in  terms  of  expected  delays.  Even  when  the  expected  delays  of  the  local  users 
approach  those  that  correspond  to  the  throughputs  of  the  LSRAAs,  the  expected  delays  of  the 
marginal  users  remain  low,  never  exceeding  ten  slots,  for  all  the  examined  cases.  The  delay 
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advantage  of  the  marginal  users,  as  compared  to  the  local  users,  increases  monotonically,  as  the 
rate  of  their  traffic  increases. 


V.  COMMENTS  AND  CONCLUSIONS 

In  this  paper,  we  studied  a  two-cluster  interconnected  system.  Each  cluster  deploys  the 
limited  sensing  random  access  algorithm  in  [4],  and  the  interconnection  is  due  to  marginal  users, 
who  dynamically  select  one  of  the  two  algorithms  for  their  transmissions.  We  performed 
rigorous  analysis  of  the  overall  system  stability,  and  the  expected  per  packet  delays.  Our  stabil¬ 
ity  analysis  generalizes  easily  to  multi-cluster  systems,  with  marginal  users  who  may  select 
either  one  of  the  local  algorithms  for  their  packet  transmissions.  For  M  clusters,  the  ergodicity 
conditions  are  as  those  in  Theorem  1,  only  that  then  there  are  M  inequalities  in  (17)  with  indices 
varying  from  1  to  M.  Then,  when  a  marginal  user  with  a  packet  to  transmit  observes  simultane¬ 
ous  beginnings  of  K  CRIs,  (from  the  algorithms  in  K  clusters,  where  2<K<M),  and  is  within  the 
examined  intervals  of  all  of  them,  he  selects  each  one  of  those  CRIs  with  probability  1/K. 

The  interconnection  policy  adopted  is  dynamic,  and  requires  no  a  priori  knowledge  of  the 
traffic  populations  and  characteristics,  and  of  the  states  of  the  involved  subsystems.  It  only 
requires  knowledge  of  the  algorithmic  rules,  and  monitoring  of  feedbacks  from  the  time  a  packet 
is  generated  to  the  time  that  it  is  successfully  transmitted.  In  addition,  the  adopted  interconnec¬ 
tion  policy  presents  a  significant  delay  advantage  to  the  marginal  users.  In  all  cases,  it  maintains 
the  value  of  the  expected  per  marginal  packet  delay  below  ten  slots,  even  when  the  expected  per 
local  packet  delays  approach  their  limit  values  within  the  stability  region  of  the  system.  This 
delay  advantage  to  the  marginal  users  may  be  of  high  importance,  when  they  transmit  high  prior¬ 
ity  data,  and  when  dynamic  cluster  reconfigurations  may  result  in  temporary  isolation  of  the 
marginal  users  if  the  transmission  of  their  data  is  delayed. 


y.*' 


Figure  7 


Expected  Per  Packet  Delays 


VI.  APPENDIX 


VI.  1.  Fundamental  Results 

We  start  by  proving  some  fundamental  algorithmic  properties  of  the  system,  which  will  be 
used  for  the  proof  of  the  theorem  and  the  lemma  in  the  main  part  of  the  paper. 

Let  us  define  the  following  quantities: 

The  number  of  CRIs  generated  by  the  algorithm  in  cluster  i, 
between  the  time  instants  T(ns)  and  T^,  where  i=l,2.  The  First 
such  CRI  starts  at  T^s\  and  the  last  such  CRI  ends  at  T„+i . 

The  length  of  the  arrival  interval  from  subsystem  k,  k=l,2,3, 
resolved  by  the  j-th  CRI  generated  by  the  algorithm  in  cluster  i, 
i=l,2,  between  the  time  points  T^s)  and  T$i.  We  note  that 
l<5jl)(k)<A;  Vi.j,  k. 

The  length  of  a  CRI,  generated  by  either  one  of  the  two  local 
LSRAAs,  which  resolves  a  length  interval  containing  arrivals 
from  an  intensity  X  Poisson  process,  and  a  length  82  interval  con¬ 
taining  arrivals  from  an  intensity  p  Poisson  process. 

L®5iE{/Oi5y>(i)  +  -y8f)(3))ld}  ;  i=l,2 
L( *  E{/(XiS<i)(i)  +  X3  5^(3))  Id)  ;  i=l,2  ,  j>2  (A.l) 

Ej|~(k)  ^  l|  j  -  E{6P(k)  Id)  ;  k=i,3,  i=l,2,  j>l 

From  the  characteristics  of  each  local  LSRAA,  (see  [4]  and  section  IV  of  this  paper),  we 
obtain: 

E{/(Xx  +  |ly)lx,y}  >  E{/(Xx+|iy)  I  x,y]  -  x> 

>  E  { / (Xx+py )  I  x,  y )  [  1  -X" 1  sup  ■- *-  -  ]  = 

x  E{/(x)} 

=  E{/(Xx+py)  lx,y][l-0.43X_1] ;  for  X>0  (A. 2) 

From  (A.l)  and  (A. 2),  we  thus  obtain: 

(X;  +  X3)  L(l]-  >  X,  E®5  +  X3  E®j  (3)  > 

-  [Xi  +  X3  -  1.29]  Lj~  ;  i=l,2 

(X,  +  X3)  >  Xi  E^  +  X3  E^  (3)  >  (A. 3) 

>[Xi+X3  -0.86]  L(j  ;  i=l,2,  j>2 


Nj: 

8j®(k): 

/  (X53  +  p.82)^ 


In  addition. 


AV(d)  =  £  E-i  (1)  +  ^2  I  E]j  (2)  + 

l  j=i  i=i 


(1)  (2)  - 
+  X3  LE^(3)  +  \3  I  EJ  j  (3)  Id 
j=i  j=i 


From  (A.3)  and  (A.4),  we  thus  conclude,  for  £  Xj  <  1: 

M 

—  n)  — 

IAV(d)  i  <max(0.43,  X!+A.3)E{2  L.-  Id}  + 


'  1  (2)  - 

+  max  (0.43,  X2+X3)  E{£L.-  Id  }  < 

j=i  J> 

N1  (1)  _  n2  (2)  - 
<E(SL ld}+E(IL  Id} 


(A-4) 


;  where. 


E{T<&  -T<s>  Id*^}  =E{£  l!j  Id}  =  E{£  lJ2-}  Id} 


We  note,  that  for  d^<*>,  Sj’^k)  =  A,  V  1,  j,  k,  and  we  then  obtain  directly  from  (A.l)  and 

0: 

lim  AV(d)  = 


=  (X1+X3)jE{/(A[X1  +  -])}-  A 


-i\l 


+  |^E{/(A[X! +X3])} -Aj  Jim  E{Nj  I d}— 1 


+  (X2  +  ^-3) j  E{/(A[X2  +  —  ])r~A 


E(/(A[X2  +  X3})}-A  Jim  E{N2  i  d}  -1 

-I  d— *»> 


;  where, 


E{/(A[Xj  +  -f])}  +  E{/(A[X,  +  X3])}[Jim  E(N]  Id}-1]  = 

2  d-*« 


=  E{/(A[?.2  +  -H)}  +  E{/(A[X2  +  X3])}[Jim  E{N2  I d}-l] 

“  d— >«« 


We  now  state  and  prove  a  lemma. 


Lemma  A 


For  all  given  Xj,  i=l,2,3: 

'  NS  ...  ]  Nj  ...  ] 

P  XL  J<xld  <P  JLj<xld'  ;  for  all  x,  for  all  d  >  d',  for  i=1.2 

j=i  J'  J  Lj=i  J> 


Thus, 


V  L(l;  Idl  >e|  TL(7ld'i  ,  foralld>d',  fori=l,2 

jS  J,d  1  J,d  J 


(A.  10) 


I  AV(d) I  <  2 Jim  E<  ^L'^ld 
U=1  J’ 


(A.ll) 


Proof 

Letd=  { d^ }  2<j<3 .  Given  d(1)  and  d(3),  consider  k  CRIs  generated  by  the  LSRA  A  in  cluster 
1.  Let  l}  denote  the  length  of  the  j-th  such  CRI,  and  let  those  lengths  be  temporarily  fixed.  Let 

N(X  /,)  denote  the  number  of  CRIs  generated  by  the  LSRAA  in  cluster  2,  within  X  h  slots. 

i=l  i=l 

Then.  since  P(/(Xd1)<x)>P(/(Xd2)^x),  if  di>di,  we  easily  conclude  that 

1c  k 

P(N(X  ld(2)=w,d(j))5J?(N(X  /i)lx  ld(2^=y,d(3)),  for  y  <  w.  Since  the  latter  inequality  is 
i=l  i=l  _ 

true  for  all  k  and  {/;},  and  since  P(/;<x ld)<P(/;£x Id");  for  all  lx,  x,  andd>d/,  we  conclude 
(A.9).  (A. 10)  is  clearly  deduced  by  (A. 9),  and  (A.ll)  results  from  (A. 10),  in  conjunction  with 
(A. 5)  and  (A. 6).  □ 

Due  to  Lemma  A,  and  (A.ll)  in  particular,  to  prove  (12)  in^Lemma  1,  it  suffices  to  show 

that  the  limit  in  (A.ll)  is  bounded,  for  all  {Xj)i<j<3  such  that  X  X,<1.  We  will  perform  the 

j=i 

latter  study  later.  At  this  point,  we  will  state  and  prove  a  proposition,  which  is  related  to  part  (ii) 
of  Lemma  1. 


Proposition  A 

^  A  3 

Given  fX, } i <j^3 ,  such  that  X_  X  <  1,  given  some  bounded  natural  number  /,  such  that 

j=i 

X/A>1,  let  us  define, 

H2^  {d^[d^]j<jS3  :  X  *-j  4(J)  <  X/A} 

j=i 


(A. 12) 


(A.13) 


Pd.U  =  ZPd.e 

e:  V(e)<V(d>-WA 

;  where  pj}^  is  as  in  (5).  Consider  F^y  (z)  in  (6).  Then, 

Fd.v  (z)  >  —  lA  —  P£u  X[V(d)-V(e)]p5.« 

e:  V(d)-V(e)>XJA 
;  for  all  z  in[0, 1),  and  all  d  in  C-H2 


(A.  14) 


Proof 


From  (6),  we  obtain: 


—  [  _XiA  1  _X_  1 

Fd.v(z)  >  (1-z)-1  zV(d)  \  -[l-PS.x./l'^r — “ '  T  + 

I  Z-1  z-1 


+  P iu-  £zv^d>  Pd,e 

Z  L  e:  V(e)-V(d)  <  -XJA 


(A.15) 


Noting  that  (1-z)-1  (l-zx)<l,  for  all  z<l  and  X<1,  and  that  (z-1)  ^“-l)  <  a,  for  all  z<l  and 
a>l,  we  obtain  from  (A.15),  the  inequality  in  (A.14).  □ 

Let  us  now  state  and  prove  another  lemma. 

Lemma  B 

A  3 

Given  {^j  )isj<3,  such  that  X  “  ]T  Aj<l,  then,  for  any  finite  natural  number  /,  and  the  result- 

j=i 

ing  then  finite  subset  H2  as  in  (A.  12),  we  have: 

P5.U  Z[V(d)-V(e)]  pi5  < 
e:  V(d)-V(e)>XM 

<  Urn  P(T&  -  T„s)  >  Id)  Z  x  P <T$i  " =  x I d) 
d-~o  A-l  * 


;  for  all  z  in  [0, 1),  and  all  d  in  C 


(A.  16) 


Proof 

Consider  the  variables  N,,  i= 1 ,2.  defined  in  the  beginning  of  this  appendix.  Then,  we 
obtain: 

P(V(d)-V(e)>LxAld)<p[  (Xi+hWi  A  +  (X2+/-ON2A-^[T^1-T[l'i)]>X.xAld] 


p\ A[N1+N'2I-Ax-[TS1-T^)]  >  Old  I  = 


*.  /.  \  A  %  /.  % 


(A.  17) 


»nrrwnrir» 


=  P 
;  where. 


T<&  -T«<A[N,+N2]-Axld 


Tn+!  -  T^s)  >  Nj  +  N2 


(A.  18) 


From  (A.  17)  and  (A.  18)  we  thus  obtain: 

P3Xx  a  P(V(d)  -  V(e)>XxA  ld)<P(N!+N2>x-^-  Id)  < 

<P(T&-T<s)>x^-ld) 

Since  (A.  19)  is  true  for  all  x,  it  also  implies, 

X[V(d)-V(e)]  p5.i  <  I  x  P  P(T«1-T<*)=xld) 


e:  V(d)-V(e)>WA 


x>l 


(A.19) 


(A.20) 


A-l 


Setting  x=  /  in  (A.19)  and  (A.20)  and  multiplying,  gives  the  first  part  of  (A.  16).  Now,  due 
to  Lemma  A  and  its  proof,  we  also  obtain: 

P(T£i  - T<s)  <  X  Id)  < P(T&  - T<s)  <  x I ?) ; 
for  d  >  d' 

which  gives  the  last  pan  of  (A.  16).  □ 

We  now  express  a  proposition,  which  ties  the  results  in  lemmata  A  and  B  with  Lemma  1,  in 
the  main  pan  of  the  paper. 


Proposition  B 


N, 


(i)  If  lim  E-j  X  L°j  ldl<°°,  for  all  {Xjh<j<3  such  that  X  ^j<L  then  (12)  in  Lemma  1  is 

d-~  J  j=l 

satisfied. 

(ii)  If  there  exist  finite  integer  l0  and  finite  positive  constant  B",  such  that , 

2 

P(Tn+i— TnS)>/-~-  I  d)  X  *  P  (Tn+fTn  )=x  ld)<B',  V/>/0,  Y{\} }  lsjs3  :XXj<l 

A-i  A  j=1 


lim 


x>i  ■ 


A-l 


then,  (13)  in  Lemma  1  is  satisfied  with  B=/0A+B'  and  H2  as  in  (A.  12),  with  /=/0.  □ 

The  proposition  is  a  direct  result  of  the  statements  in  Lemmata  A.  and  B.  Let  us  now  con¬ 
sider  the  variables  Ni({Xj},  A),  i=l ,2,  defined  in  Section  IV.l  of  the  paper,  and  the  generalized 
drifts  A(i)V(dA,  (Xj}),  i=l,2,  and  AV(dA,  {Xj}),  in  (14)  and  (15)  respectively.  Towards  the 
proof  of  Theorem  1,  in  Section  IV.l,  we  then  state  and  prove  a  lemma. 

Lemma  C 

Let  the  variables  N,({X,  },A),  i=l,2,  be  such  that: 


a.5 


There  exist,  finite  natural  number  n<j,  p:  0<p<l,  and  finite  positive  constant  c,  such  that: 

P(Ni({Xj},  A)=n)  <cpn,  Vn>n0,  V{Xj}:  £  Xj  <1,  i=l,2  (A.22) 

j=l 

Then, 

(i)  If  the  conditions  (17)  in  Theorem  1  are  satisfied  for  some  E>0,  then,  there  exists  a  finite 
subset  of  C,  and  some  £'>0,  so  that  (9)  is  satisfied.  That  is,  the  Markov  chain  {Sn},£0 
is  then  ergodic. 

(ii)  If  at  least  one  of  the  two-generalized  drifts  A(l)V(dA,  {X-j }),  i=l,2,  is  nonnegative,  then 
there  exists  a  proper  subset  H  of  C,  so  that  (11)  is  then  satisfied.  That  is,  the  Markov 
Chain  {Sn  1,^0  is  then  nonergodic. 

Proof 

(i)  Let  (A.22)  be  true,  and  let  Y  X:  <1.  Define  then,  sup  I  Ex{/IA}-A  I.  Given  £  in 

pi  ~  X<1  , 

(17),  select  l0  such  that  l0  >  n0  andac  £  pn  =  -^-  =  (l-p)“1  acp/a+1;  that  is,  select 

n>/„+l  2 

lQ  =  max(n0+l,  [/np]_1/n(£[l-p][2ac]-1  -1).  Select  then,  Hj  ^  {d^  [d^]i<j<3  : 

d(1)</0,  d®</0,  d(3)<2/0}.  Then,  due  to  (17),  and  in  conjunction  with  (A.9)  in  Lemma  A 
—  £  — 

we  conclude:  AV(d)<-y,  for  all  din  C-Hj.  The  proof  of  part  (i)  in  the  Lemma  is  now 
complete. 

(ii)  Let  (A.22)  be  true,  and  let  £  Xj  -el  .  Let  AwV(dA,  {Xj  })>0,  for  either  i=l  or  i=2.  If  then, 

M  - 

there  existed  Hj  subset  of  C,  defined  as  in  the  proof  of  part  (i),  such  that  AV(d)<-£,  for  all 

C-Hi  and  some  e>0,_then  due  to  (A.22),  the  constant  lQ  in  Hi  could  increase  (remaining 
finite),  to  give  A(l)V(dA,  {Xj  })<0.  Thus,  part  (ii)  of  the  Lemma  is  true.  □ 

As  it  is  apparent  from  the  statements  of  Proposition  B  and  Lemma  C,  to  complete  the 
proofs  of  Lemma  1  and  Theorem  1,  in  Section  IV.  1,  we  need  to  study  the  statistical  behavior  of 
the  variables  involved  in  the  description  of  the  time  interval  [T^*  T$i],  when  the  lags  at  T„s) 
are  sufficiently  long,  so  that  all  CRIs  in  [T^s) ,  T^  ]  resolve  arrival  intervals  of  length  A,  from  all 
the  three  subsystems.  We  perform  such  studies  below. 


VI. 2.  Studies  for  Sufficiently  Long  Initial  Lags 

Our  objective  here  is  to  study  the  behavior  of  the  variables  N;({Xj},A),  i=l,2, 
A(l)V(dA,  (Xj }),  i=l,2,  and  AV(dA,  { Xj }),  defined  in  Section  IV.  1.  For  tractability,  we  introduce 
simpler  notation.  Give  {X}  J  i<j<3,  we  first  define: 

Xj  t  (Xj+X3)A,  i=l,2 


-A  -J* 


y  ^  (X3A)I2 


(A.23) 


MiZ^+X3,Pl,2 

For  the  rest  of  this  subsection,  we  assume  that  every  CRI,  within  the  interval  [T^*  t£2i], 
resolves  arrivals  in  a  window  of  length  A,  from  all  the  relative  subsystems;  that  is,  a  CRI  gen¬ 
erated  by  the  algorithm  in  cluster  i,  resolves  a  A-length  window  from  subsystem  i,  and  a  A-length 
window  from  subsystem  3.  We  then  define: 

i=l,2,  ;  Pi(e,m):  The  probability  that  T^s)+e  is  a  collision  resolution 

point  for  the  LSRAA  in  cluster  i,  and  that  it  is  the 
m-th  such  point  We  note  that  m=l  refers  to  the  end 
of  the  first  CRI,  which  starts  at  T^,s). 

i=l,2  ;  P(l)(e,m):  The  probability  that  T$j  -  Tj,s)  =  e,  and  there  are  m 

CRIs  for  the  algorithm  in  cluster  i,  in  [T$f\  T^iJ. 

i=l,2  ;  Pw(m):  The  probability  that  there  are  m  CRIs  for  the  algo¬ 

rithm  in  cluster  i,  in  [T^s),  T$j]. 

/* :  The  length  of  a  CRI,  when  it  resolves  an  interval  of 

length  d  containing  arrivals  from  an  intensity  A.  Pois¬ 
son  process,  where  A.  and  d  are  such  that,  Ad=x. 

0<«<k  ;  /n-k_n:  The  number  of  slots  needed  by  either  one  of  the  two 

LSRAAs  in  the  system  to  transmit  k  packets,  when  n 
of  them  have  counter  values  equal  to  one,  and  the 
remaining  k-n  packets  have  counter  values  equal  to 
two. 

0</<e-li-  ;  Pi(e,m;/,n):  The  probability  that  T|,s)+e  is  a  collision  resolution 

i=1-2  J  point  for  the  LSRAA  in  cluster  i,  that  it  is  the  m-th 

such  point,  that  T^\  -  T„ : ’>e,  that  the  last  before 
T^s)+e  collision  resolution  point  for  the  LSRAA  in 
the  other  cluster  occurs  at  T^s)+  /,  and  that  it  is  the 
n-th  such  point. 

i,  =  1,2  ;  Pitic(e,m;e,n):  - -rts)  _-r(s)^  _ 


^  The  probability  that  =  T^+e,  and  that  this  point 

>c  =  |L  ifi=2  is  the  m-th  collision  resolution  point  for  the  algo- 

(_2,  if  i=l  nthm  in  cluster  i,  and  it  is  the  n-th  such  point  for  the 

algorithm  in  cluster  ic. 

To  stan  with,  the  deployed  LSRAAs  induce  the  following  recursions,  regarding  the  lengths 
fn.k-n  (see  [4]): 

P(A.k-t  =  m)  =  Pf/k-t.o  =  m-1) 

P(/0,k=m)  =  P(4.0  =  m-1) 

P(/,,o  =  1)  =  P(/0,o  =  1)  =  P(/0,i  =  2)  =  P(/u  =  2)  =  1 


I' 

<* 

.* 

>, 

> 

i 


c 


fr 


I 


k>2 
m>k+l 

k>n 

n>3 

m>2k-l 


;  P(/zk-2=m)  =  2~2  ]  P(/k,0=m-2)+P(/2.k-2=m-l)+2P(/k-i,o)  =  m-2) 


;  P(/n,k-n=m)=2"ni  P(^o=m-2)+P(/iuk-n=ni-U+nP(/k-i.o=m-2)+ 


n-l  r  „ 

*5l: 


P(/i.k-i=m-l)| 

In  addition,  for  J_J  denoting  integer  part,  we  have: 


(A.24) 


— r  X 


P(/X=m)=  I  e~x—  P(/k.o=m) 

We  also  have  the  following  recursions  and  expressions: 

P(/X._y=e)  ;  if  m=l 


(A.25) 


'eim  j  ’  p'^e,m) =' 


£  Pj(k,m-l)P(/x.  =  e-k) ;  if  m>2 
fc=m-l 


(A.26) 


For, 


we  also  have: 


.  A 


1  ,  if  i=2  A  f  1,  if  x>0 

2  ,  if  i=l  U(x)=j  0,  ifx<0 


K 

f. 

V, 

i=l 

u  . 

4 

i 

m>2 

*r. 

V" 

n>l 

1 

e>m 

✓ 

[>n 

•p. 

i=l,2 

U* 

i 

i=l,2,  m>l,  e>m;  Pi(e,m;0,0)  =  Pi(e,m)P(/Xie_y>e) 


m>2,  e>m"l 

o  !»  ;  Pi  ic(e,m;e,l)=Pi(e,m)P(/X]c_y=e) 
i  i,—  I  - 


(A.27) 


(A.28) 


(A.30) 


i  min(e-M,e-m+t)  .  _ ,  ,, 

[■  ;  Pi(e,m;/,n)  =  £  Pi(e-k,m—  l;/,n)P(/Xj=k)P(/Xjt  >e-[)P^  (/Xjt>e-k-0  + 

k=l 


+U(7-m)  *5?"  Pic(/,n;e-k.m-l)P(/Xj=k)P(/x.  >e-/)P_1(/x  >/+k-e)(A.31) 
k=«-/+i 


a.8 


% 


m>4 


n^2^i.ie(e,m;e,n)=  2 


e-l  min(e-M,e-t-l-m) 


Pi  (e-k,  m-1 n- 1  )P(/Xt  =k). 


J 


f=n+l  lc=l 


.P(/x.=e-/)P‘1(/xi,>e-k-/) 


e-l  min(e-/-l.e+l-n) 

+  2  2  (e-k,  n- 1 m- 1  )P(/X.  =k)P(/Xj  =e-l)p- 1  (/x,  >e-k-/)  (A.32) 

i=m-l  k=l 


P(l)(e,m)  =  2  Pu  (e,m;e,n)  ;  i=l,2 


(A.33) 


P®(m)  =  2  P{i)(e,m) ;  i=l,2 

e2m 


(A.34) 


Referring  to  the  expressions  in  (14)  in  Section  IV.  1  of  the  paper,  and  from  (A.33),  we  also 
obtain: 


i=l,2;  A®V(dA,  {Aj })  =  2  2 P(i)(e,m)[e-mA] 

m>l  e^m 


=  2  IP(i)(e,m)(e-m— ] 
m21 


(A.35) 


Subject  to  the  satisfication  of  Lemma  1  and  Theorem  1,  in  Section  IV.  1,  the  h<j<3 
values  for  which  the  Markov  Chain  {S„}„20  is  ergodic,  are  found  as  those  that  make  the  expres¬ 
sions  in  (A.35)  negative. 

Via  some  tedious  manipulation  of  the  expressions  found  in  this  section,  which  we  do  not 
include  here,  we  found  the  result  stated  in  the  proposition  below. 


Proposition  C 

Given  e>0,  there  exist  positive  integer  m,,  and  0<|I<1,  such  that: 
i=l,2;  2  +  H  2  P®(0 


fem+l 


l >m 


Vm>mo  ,  V-  {Xj  }lsjs3  :  2  Xj  <  1  (A.36)  □ 

We  are  now  ready  to  complete  the  proofs  (flemma  1  and  Theorem  1,  in  Section  IV.  1.  We 
do  this,  in  Section  VI.3  below. 


VI.3  The  Proofs  of  Lemma  1  and  Theorem  1 

From  Proposition  C  and  expression  (A.36),  we  conclude  that  there  exist  constants  n0,  p, 
and  c,  such  that  (A.22)  in  Lemma  C  is  satisfied.  Thus,  the  statements  in  the  latter  Lemma  hold. 

In  addition,  due  to  (A.22),  (A.21)  in  Proposition  B  clearly  holds  as  well  then,  and  also 
N\  ^  _  3 

JimE{2  L.  -  ld}o>,  for  all  {Xj)i<jS3  such  that  2  ^j<l-  Thus,  the  statements  in  Proposition  B 

d— *<»»  j=[  ■*'  j=r  1 

and  in  Lemma  C  all  hold,  and  so  do  Lemma  1  and  Theorem  1 . 


c 


0 


4 


% 
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VI.4  Recursions  and  Systems  for  Delay  Analysis 


Here,  we  present  the  quantities  needed  for  the  computation  of  the  expected  per  packet 
delays,  and  their  recursions.  Let  us  define,  for  4  as  in  (A.27): 


(di.Ss;  {/,  dje,  d3}): 
;  i=l,2 


hi(dj,  53;{Uie,d3}): 

;  t=l,2 


^(di,  53;{/,dit,d3}): 

is  1,2,  j=  1,2,3 


(di,S3;{/,dic,d3}): 

;i=  1.2,  j=  1.2,3 


0p)(di,53;{/,dit,d3}): 

;  i  =  1,2  ,  j  =  i.3 


wp)  (di,  83;  {/,  d^,  d3 }) : 
i=  1,2  ,  j=  1,3 


^n,k-n  •’ 
;  nSO ,  k>n 


During  the  process  of  the  algorithmic  system,  the  event  of 
being  at  some  time  instant  t  when  the  first  slot  of  a  CRI  for 
the  LSRAA  in  cluster  i  begins,  that  at  t  the  lag  for  subsys¬ 
tem  i  is  d;  and  the  lag  for  subsystem  3  is  83,  that  the  last 
before  t  CRI  for  subsystem  ic  starts  at  t and  that  at  t -/  the 
lag  for  subsystem  ic  is  d^  and  the  lag  for  subsystem  3  is  d3 . 

Given  the  event  (d{,  S3;(/,  d^,  d3 }),  the  length  of  the  time 
interval  until  the  first  after  that  occurence  of  some  T^s)  with 
Sn  equal  to  the  unity  vector. 

Given  the  event  (£^,83;  {/,  d^,  d3 }),  the  expected  cumula¬ 
tive  delay  of  all  the  packets  from  subsystem  j  that  are 
transmitted  from  the  instant  when  the  above  event  occurs, 
to  the  instant  when  the  first  after  that  T„s)  with  Sn  equal  to 
the  unity  vector  occurs. 

Given  the  event  (d,,  S3;  {/,  d^,  d3 }),  the  expected  cumula¬ 
tive  delay  that  those  packets  from  subsystem  j,  which  are 
waiting  at  the  point  of  the  event,  have  already  experienced. 

Given  the  event  (di,  83;  {/,  dic,  d3}),  the  expected  cumula¬ 
tive  waiting  time  of  those  packets  from  subsystem  j  which 
are  transmitted  during  the  CRI  of  the  LSRAA  in  cluster  i 
that  starts  with  the  event. 

Given  the  event  (d,,  S3;  {/,  d^.,  d3 }),  the  expected  cumula¬ 
tive  transmission  time  of  these  packets  from  subsystem  j 
which  are  transmitted  during  the  CRI  of  the  LSRAA  in 
cluster  i  that  starts  with  the  event. 

Given  that  during  some  CRI  from  any  of  the  two  LSRAAs, 
n  packets  have  counter  values  equal  to  1  and  k-n  packets 
hae  counter  values  equal  to  2,  the  expected  cumulative 
delay  of  all  the  packets  transmitted  from  the  point  of  the  (n, 
k-n)  occurence  to  the  end  of  the  CRI. 


To  start  with,  the  operations  of  each  of  the  two  LSRAAs  give  rise  to  the  following  recur¬ 
sions: 


Zq.o  -  0 .  Z1.0  -  L  Zo.k  =k  -f-Z^o  ;  k>l 
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1 


J* 


(A.37) 
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R 
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r* 

> 
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Zu_i  =  k  +  Zo,fi  =  2k- 1  +  Zfc_i.o  ;  k£2 

z^-k+r-ff?]  Zi.k.ijnaz 

i=0  U  J 


Then,  defining 


X3  ;  if 1 
V3  =1  X3 


Cj  = 


-j- ;  if /=  0 


JXj  min  (dj,  A);ifj=i 
^  v3min(83,  A);j=3 


(A.38) 


A.39) 


we  have: 


i  =  1,2 ;  wp>  (di,  83  ;  {/,  di, ,  d3 })  =  w®  (di,  83)  = 

^j[Xj  min(di,  A)  +  v3  min(53,  A)]-1,  exp{-Xj  min(di,  A)-v3  min^,  A)} 
[Xj  min(di.  A)  +  v3  min(83,  A)]k 


s 

kao 


k! 


-to 


(A.40) 


In  addition, 

i— 1,2 ;  0®  (1,  83;  {/,  d,,.  d3})  =  Ap  (1,  83  ;  {/,  di,,  d3)=  1.5^  ;  ¥83,  {/, d^,  d3} 

i=  l,2;0P)(di>  1;  {/,  dj,.  d3})=A|3>(di,  1;  U,\,  d3}>=  1.5X3;  V  l>  l  (A.41) 
0P  (di,  83;  {0,  di,,d3 })  =  A®  (di,  83;  {0,  d,.,  d3 })  =0,  for  83*  d3 
AP  (di,  d3;  (0,  d2  ,  d3})  =  A$>  (d2,  d3;  {0,  dj,  d3});j=l,2,3 
i=  1,2;  AP>  (di,T  ;  {O.d^  ,  1})  =  2  ©P>  (di,  1 ;  {0,  d^,  1})=  1.5X3  ;  ¥di,  d^ 

i=l,2 ;  AP  (di,  83  ;  {/,  d^,  d3J)  *0,®  (di,  83  ;  {/,  dj,,  d3J) ;  ¥di  £  A,  ¥83,  {/,  dj.,  d3} 


i=l,2;  9p5  (di,  83 ;  (/,  d^,  d3 })  =•< 


AP>  (di,  83  ;  {/,  dj,,  d3}) ;  V83<A,  if  fitl 
2-1  AP>  (di,  83 ;  {/,  di.,  d3 }) ;  V  83^,  if  h 0 


i=l,2;  A-1*5  (di,  83;  {/,  i,d3})=-| 


1.5Xi,  ;  if /  =  0 

/(4  +  l)Xi,;if/>1,Vdi,53,d3 
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For  V3  as  in  (A.39)  and  Cj  as  in  (A.40): 

^3}  •.  a,®  (di,  53 ;  {/,  di,.  d3 })  -  0©  (di,  53 ;  {/,  dj.,  d3 })  = 

A P (dj+m-Xf 1  £j , 83  +m-min(53 , A); { /+m, d^ , d3  })-m[Xj  dj— 

-Xjm(~t-l);w.p.  P(/^+Vl=m)P(/^+v,  >m+0 

A  P  (dj+m-Xf1  Ci  ,53+m-min(83 ,  A);  { 0,  d^  +m+/-X“1  ,  83+m-min(83 ,  A) } ) 

-m[Xj  dj-C j  ]— Xj  m(  y+1  );  w.p.  P(/<r+Vj  =m)P(^  +v,  =m+0 
A[ci)(die+m-Xi;1Cie,  d3+m-min(d3,A);{m-/,di,S3 ))-(m-/)[Xj  Cj-dj] 
-Xj(m-/)[^— +1]»  for  m>/;w.p.  P(/^+Vj=m)P(/^.+V3>m-/) 

hi  (di,  d3  ;  {0,  d2,  d3 })  =  h2(d2,  d3;  {0,  di,  d3 )) 

For  d;  <  A,  d^  <  A,  d3  <  A ;  h^d,,  d3  ;{0,  d^,  d3 })  =  1 ;  w.p.  P (/^  =  l)P(/j;jc+^3=l) 
i=l»2;  hi(di,83;{/,die,d3})  = 

m  +  hiCdi+ra-X^Ci,  S3+m-min(83,A);{/+m,  die,  d3 }) 

;  w.p.  P  (/^i+V3=m)  P(/^+v3  >  m+0 

ra  +  hiCdi+m-Xf1  Ci.  83  +  m-min(83,A);  {0,  dit+nM-/-X^1Cic,S3+m-min(53,A)}) 
;  w.p.  P(/(r+v3=m)P(/^+vj=na+0 

m  +  (d^+m-X^1  Ci^  .  ^3  +  m-min(d3 , A); {m-/,di,  83 }),  for  m  >  /, 

;  w.p.  P(/^+V}=m)P(/);j+Vj  >  m-0 

We  note  that  the  expected  value  H  in  ( 1 8)  is  such  that,  H  =  hi  ( 1 ,  1 ;  {0,1,1}) 

=  h2(  1,1;  {0,1,1  }). 

Defining, 

0P  (d; , 83 )  =  I  0p)(di,53;{/,dit,d3})P(di,83;{/,dic,d3}  ^,83) 
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we  also  obtain: 


Cl 


wP  (l,l)  +  eP(l,l;  {0,1,1}),  forj=l,2 
;  w.p.  P(/x,+x,  i2=l)P(^+X,  12  =  1) 
i  =  1,2  vP  (1, 1;  {0, 1, 1 })  =■  wf>  (1,1)  +  w£>  (1,1)  + 

+  0\3>  (1,1;  {O,l,l})  +  0£>  (1,1; {0,1,1}),  forj=3 

;  w.p.  P(/x,  +  x,  i2  =  mix,  +  X,  12  =  1) 


(A.45) 


vP  (di,  83;  {0,  die,  d3})  =  0,  if  83  *d3,  i=l,2,  j=l,2,3 


'wP(di,d3)  +  eP(di,d3;{0,die,d3}) 

,  for  j=l,2;  w.p.  P(/xidi+d1X,i2=l)P(/xitdie+d3^i2=l) 


vP(di,  d3;  {0,  d^,  d3})=«  w(3)  (dj,  d3)  +  wf  (d^,  d3) 

:  for  d,  5  A.  4,  5  A.  d,  <  A  +  ep)(dij  {a  dui  d3 })  +  0(3)  (d.c>  d3 ;  {O.dj,  d3 }), 


[,  for  j=3;  w.p.  P^+^x,  i2=l)P(/^dit+<i,X,  12=1) 
vP  (di,  d3;  {0,  d2,  d3}=yP(d2,  d3;  {0,  di,  d3}) ;  j=l,2,3 


(A.46) 


wP(di,53)  +  eP(di,83;{/,die,d3}) 
vP ^(di . 83 ; { /, d^ , d3 } )  =  +vp(d.+m_xri;i,  53+m-min(53,A) 
j=l,2  ;{/+m,dic,d3}) 


;  w.p.  P(/^+Vj=m)  P(/^+v,  >  m+/) 
wP(di,53)  +  0P(di,53;{/,d^,d3}) 

+  ¥P(dj+m-Xj'1  ^,53+ra-min(83,A); 
{Q.d^+m+Z-X^1^,  53+m-min(53,A)}) 
;  w.p.  P(/?i+V3=m)  P(/^+Vj=m+/) 


(A. 47  a) 


d3+m-tnin(d3.  A); 

yP  (di,53;{/,dic,d3})=v 

j=i, 2  • for  m>/;  W-P-  p(^  +v3  =m)P(/^i+v,  >m-/) 


(A.47b) 


vP)(di,53;{/,die.d3}) 


vP)(di,63;{/,di,,d3}) 


wP)(di,53)  +  e(3)(di,53;  {/,di.,d3}) 

+  yP^di+m-Xr1  ^,53+m-min(53,A);{/+m,dic  ,d3 }) 
;  w.p.  P(/^+Vj=m)  PC/^+v,  >  m+/) 

wP)(di,S3)  +  eP)(di,83;[/,dii43D 
+wp>(die,d3)  +  e(c3>(dis,d3)  + 

\yP  (dj +m-Xp  ^,§3  +m-min(S3 ,  A; 

{ O.d^ +m+/-X^  Ci, ,  63  +m-min(S3 ,  A) } ) 

;  w.p.  P(/£+v,=m)  P(/^+Vj=m+/) 

w 

wP>(dii,<b)  +  ef)(di,,d3) 

+  Vif(die+m-A.jp  £ic,d3+m-min(d3,A);{m-/,di,53 }), 
for  m>/; 

;w.p.  P(/?.  +Vj=m)P(/?,+Vj  >  m-/) 


(A.48a) 


(A.48b) 


We  note  that  the  expected  value 

W<j)  =  v|/P(U;{0,U})  =  VP(U;[0,U}),j=l,2,3. 


The  relationships  included  in  this  section  induce  infinite  dimensionality  linear  systems.  We 
used  the  methodology  in  [6]  to  derive  upper  and  lower  bounds  on  those  systems,  and  subse¬ 
quently  upper  and  lower  bounds  on  the  delays  D^\  j=l,2,3. 


a.  14 


REFERENCES 


[1]  A.  Ephremides,  J.  E.  Wieselthier,  and  D.  T.  Baker,  "A  Design  Concept  for  Reliable  Mobile 
Radio  Networks  with  Frequency  Hopping  Signaling,"  Proceedings  of  the  IEEE,  vol.  75,  pp. 
56-73, Jan.  1987. 

[2]  B.  Eklundh,  "Channel  Utilization  and  Blocking  Probability  in  a  Cellular  Mobile  Telephone 
System  with  Directed  Retry,"  IEEE  Trans.  Comm.,  vol.  COM-34,  pp.  229-337,  April  1986. 

[3]  L.  Georgiadis  and  P.  Papantoni-Kazakos,  "A  0.487  Throughput  Limited  Sensing  Algo¬ 
rithm,"  IEEE  Trans.  Inf.  Th.,  vol.  IT-33,  pp.  233-237,  March  1987. 

[4]  M.  Paterakis  and  P.  Papantoni-Kazakos,  ”A  Simple  Window  Random  Access  Algorithm 
with  Advantageous  Properties,"  Proceedings  of  INFOCOM’88.  Also,  submitted  for  journal 
publication. 

[5]  M.  Paterakis,  L.  Georgiadis,  and  P.  Papantoni-Kazakos,  "On  the  Relation  between  the  Fin¬ 
ite  and  the  Infinite  Population  Models  for  a  Class  of  RAAs,"  IEEE  Trans.  Comm.,  vol. 
COM-35,  pp.  1239-1240,  Nov.  1987. 

[6]  L.  Georgiadis,  L.  Merakos,  and  P.  Papantoni-Kazakos,  "A  Method  for  the  Delay  Analysis 
of  Random  Multiple  Access  Algorithms  whose  Delay  Process  is  Regenerative,"  IEEE  Jour¬ 
nal  on  Selected  Areas  in  Communications,  vol.  SAC-5,  pp.  1051-1062,  July  1987. 

[7]  W.  Szpankowski,  "Stability  Conditions  for  Multidimensional  Queueing  Systems  with 
Applications,"  J.  Operations  Research,  to  appear. 

[8]  A.  G.  Pakes,  "Some  Conditions  for  Ergodicity  and  Recurrence  of  Markov  Chains,"  Oper. 
Res.,  vol.  17,  pp.  1058-1061,  1969. 

[9]  W.  Szpankowski,  "Some  Theorems  on  Instability  with  Applications  to  Multi  Access  Proto¬ 
cols,"  J.  Operations  Research,  to  appear. 

[10]  R.  L.  Tweedie,  "Criteria  for  Classifying  General  Markov  Chains,"  Adv.  Appl.  Prob.,  vol.  8, 
pp.  737-771,  1976. 

[11]  M.  Kaplan,  "Sufficient  Condition  for  Nonergodicity  of  a  Markov  Chain,"  IEEE  Trans.  Inf. 
Th.,  vol.  IT-25,  pp.  470-471,  1979. 


V^V.YVS 


■V: 


Copy  No , 
1  -  6 


DISTRIBUTION  LIST 


Director 

Naval  Research  Laboratory 
Washington,  D.C.  20375 

Attention:  Code  2627 

7  Dr.  R.  N.  Madan 

Electronics  Division,  Code  1114SE 
Office  of  Naval  Research 
800  N.  Quincy  Street 
Arlington,  VA  22217-5000 

8  Mr.  James  G.  Smith 
Office  of  Naval  Research 
Code  1241 

800  N.  Quincy  Street 
Arlington,  VA  22217-5000 

9  Professor  Mike  Athans 
MIT,  Building  35 
Cambridge,  MA  02139 

10  Professor  A.  Ephremides 
Electrical  Engineering  Dept. 
University  of  Maryland 
College  Park,  MD  20742 

11  Dr.  Dave  Castanon 
ALPHATECH,  Inc. 

2  Burlington  Executive  Center 
111  Middlesex  Turnpike 
Burlington,  MA  01803-4901 

12  Dr.  John  P.  Lehoczky 
Department  of  Statistics 
Carnegie  Mellon  University 
Schenley  Park 

Pittsburgh,  PA  15213-3890 

13  Professor  Alex  Levis 
MIT,  Building  35 
Cambridge,  MA  02139 

14  Dr.  Harold  L.  Hawkins 
Manager,  Perceptual  Sciences 
Code  1142 

Office  of  Naval  Research 
800  N.  Quincy  Street 
Arlington,  VA  22217-5000 


(Continued  on  Next  Page) 


15 


17 


29 


33 


Dr.  G.  Malecki 
Office  of  Naval  Research 
Code  1142 

800  N.  Quincy  Street 
Arlington,  VA  22217-5000 

16  Dr.  David  Kleinman 

ESE  Department,  U-157 

The  University  of  Connecticut 

STORRS,  CT  06268 

-  28  Defense  Technical  Information  Center,  S47031 

Building  5,  Cameron  Station 
Alexandria,  VA  22314 


-  30 

P. 

Papantoni -Kazakos ,  EE 

31 

M. 

Paterakis 

32 

L. 

Ming 

-  34 

E. 

H.  Pancake,  Clark  Hall 

35  SEAS  Publications  Files 


Office  of  Naval  Research  Resident 
Representat ive 
818  Connecticut  Avenue,  N.W. 

Eighth  Floor 
Washington,  D.C.  20006 

Attention:  Mr.  Michael  McCracken 

Administrative  Contracting  Officer 

Send  Cover  hotter  Only 


JO:-- 1  7  36  :  j  lb 


UNIVERSITY  OF  VIRGINIA 
School  of  Engineering  and  Applied  Science 

The  University  of  Virginia’s  School  of  Engineering  and  Applied  Science  has  an  undergraduate 
enrollment  of  approximately  1,500  students  with  a  graduate  enrollment  of  approximately  560.  There 
are  150  faculty  members,  a  majority  of  whom  conduct  research  in  addition  to  Caching. 

Research  is  a  vital  part  of  the  educational  program  and  interests  parallel  academic  specialties. 
These  range  from  the  classical  engineering  disciplines  of  Chemical,  Civil,  Electrical,  and  Mechanical 
and  Aerospace  to  newer,  more  specialized  fields  of  Biomedical  Engineering,  Systems  Engineering, 
Materials  Science,  Nuclear  Engineering  and  Engineering  Physics,  Applied  Mathematics  and  Computer 
Science  Within  these  disciplines  there  are  well  equipped  laboratories  for  conducting  highly  specialized 
research  All  departments  offer  the  doctorate;  Biomedical  and  Materials  Science  grant  only  graduate 
degrees.  In  addition,  courses  in  the  humanities  are  offered  within  the  School. 

The  University  of  Virginia  (which  includes  approximately  2,000  faculty  and  a  total  of  full-time 
student  enrollment  of  about  1 6,400) ,  also  offers  professional  degrees  under  the  schools  of  Architecture, 
Law.  Medicine,  Nursing,  Commerce,  Business  Administration,  and  Education.  In  addition,  the  College 
of  Arts  and  Sciences  houses  departments  of  Mathematics,  Physics,  Chemistry  and  others  relevant 
to  the  engineering  research  program.  The  School  of  Engineering  and  Applied  Science  is  an  integral 
part  of  this  University  community  which  provides  opportunities  for  interdisciplinary  work  in  pursuit 
of  the  basic  goals  of  education,  research,  and  public  service. 


